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Abstract 


High  performance  propellants  and  explosives  are  based  on  the  binding  of  energetic  crystalline 
solids  within  a  cured  elastomeric  polymer  matrix  in  order  to  impart  acceptable  insensitivity 
and  good  mechanical  properties  to  the  material.  Most  plastic-bonded  explosives  (PBX)  used 
are  based  on  the  inert  rubbery  hydroxyl-terminated  polybutadiene  polymer  (HTPB).  The  total 
energetic  density  of  PBXs  may  then  be  enhanced  by  replacing  this  inert  polymer  by  an 
energetic  one.  Mechanical  properties  of  the  material  should  be  remained  acceptable  with  the 
use  of  an  appropriate  plasticizer.  The  glass  transition  temperature  named  T  g  of  a  poly  mer 
blend  is  a  representative  data  and  an  useful  measurement  which  allows  to  assess  how  the 
polymer  “rubbery”  is,  and  as  well  how  the  material  sensitivity  is.  This  T  g  property  was 
predicted  by  computer  modelling,  and  was  also  measured  experimentally  for  inert  HTPB  and 
energetic  GAP  polymers. 


Resume 


On  formule  les  propergols  et  les  explosifs  de  haute  performance  en  incorporant  des  substances 
energetiques  cristallines  dans  une  matrice  de  polymere  elastomeriques  afin  de  leur  impartir 
une  insensibilite  et  des  proprietes  mecaniques  acceptables.  La  plupart  des  explosifs  plastiques 
contiennent  du  polybutadiene,  un  elastomere  inerte.  Or,  la  densite  energetique  totale  de  ces 
materiaux  explosifs  a  liant  plastique  pourrait  etre  accrue  en  substituant  le  polyazoture  de 
glycydile,  une  substance  energetique,  au  polybutadiene  a  terminaison  hydroxy  inerte. 
Toutefois,  on  devrait  employer  des  plastifiants  appropries  pour  conserver  les  proprietes 
mecaniques  du  materiau.  La  temperature  de  transition  vitreuse,  Tg,  est  un  parametre  important 
permettant  d’evaluer  les  proprietes  elastomeriques  du  materiau  et,  ainsi,  son  insensibilite. 
Nous  avons  calcule  les  Tg  du  polybutadiene  h  terminaison  hydroxy  et  du  polyazoture  de 
glycydile  par  modelisation  et  les  avons  mesure  en  laboratoire. 
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Executive  summary 

The  search  of  new  energetic  materials  more  powerful,  less  sensitive  or  more  benign  for  the 
environment  is  a  very  active  field  of  research.  Indeed,  the  army  of  tomorrow  will  have  to  face 
challenges  which  will  require  of  energetic  materials  (EM)  fulfilling  specific  constraints.  For 
example,  the  increase  in  the  range  of  a  weapon  will  ask  for  the  discovery  of  more  powerful 
EM  with  higher  performance  in  term  of  delivered  energy.  Moreover,  these  new  ingredients 
will  have  to  meet  increasingly  strict  standards  of  vulnerability  without  forgetting,  that  to 
ensure  sustainable  training  of  the  troops,  these  new  molecules  will  have  to  be  non  toxic  and 
without  effect  on  the  environment. 

Such  challenges  will  be  difficult  to  meet  since  the  requests  seem  sometimes  contradictory.  For 
instance,  we  ask  for  more  performance  as  well  as  les  sensitivity.  To  face  such  challenges  it  is 
thus  imperative  to  develop  tools  which  will  facilitate  our  task  and  will  allow  us  to  make  a 
selection  among  several  candidates.  Indeed,  the  processing  of  new  energetic  materials  is  an 
expensive  task  which  takes  time  and  moreover,  generates  an  impressive  quantity  of  waste 
which  will  have  to  be  eliminated  while  respecting  increasingly  strict  and  expensive 
environmental  standards.  Furthermore,  all  these  efforts  could  be  useless  if  the  processed 
products  do  not  meet  the  expected  characteristics. 

It  would  be  thus  interesting  to  be  able  to  evaluate  a  priori  certain  properties  of  an  EM. 

It  is  particularly  pertinent  if  one  could  predict  some  mechanical  properties  of  binder- 
plasticizer  blend.  Such  a  property  is  the  glass  transition  temperature.  In  fact,  this  property  is 
closely  related  to  the  range  of  temperature  an  energetic  composition,  such  as  a  rocket  motor  or 
a  plastic  bonded  explosive,  would  keep  proper  mechanical  properties.  Such  as  a  the  prior  to 
evaluation  of  the  glass  transition  temperature  would  avoid  spending  time  and  money  on 
useless  products  and  allow  us  to  select  only  the  most  promising  materials 

This  work  precisely  describes  the  efforts  carried  out  with  the  RDDC  Valcartier  by  using 
computational  chemistry  in  order  to  predict  the  glass  transition  temperature  for  two  very 
popular  and  used  binders:  hydroxyl  terminated  polybutadiene  and  the  polyglycidylazide. 

This  innovating  work  opens  the  door  for  the  molecular  simulation  which  will  be  a  new  tool 
for  defense  science  and  which  will  prove  surely  more  useful  in  the  future  to  meet  the  needs  of 
the  army  of  tomorrow  in  fields  as  varied  as  the  fire  power,  the  vulnerability  and  the  safety  of 
munitions. 


Abou-Rachid,  H.,  Lafleur-Lambert,  X.  Lussier,  L.-S.,  Mounir  Jaidann,  Gagnon,  N.  and  Brisson,  J. 
2008.  Effet  des  plastifiants  sur  la  temperature  de  transition  vitreuse,  Tg,  du  butadiene  &  terminaison 
hydroxy  et  du  polyazoture  de  glycidyle  —  Simulations  et  experiences.  TR  2007-615  RDDC 
Valcartier. 
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Sommaire 


La  mise  au  point  de  nouveaux  materiaux  energetiques  plus  puissants,  moins  sensibles  et 
moins  nocifs  pour  l’environnement  est  un  domaine  de  recherches  tres  actif.  En  effet,  1’armEe 
de  demain  sera  confrontee  a  des  defis  qui  exigeront  des  materiaux  energetiques  satisfaisant  a 
des  contraintes  precises.  Par  exemple,  l’allongement  de  la  portee  d’une  arme  commande  la 
decouverte  de  materiaux  Energetiques  a  la  fois  plus  puissants  et  liberant  plus  d’Energie.  En 
outre,  ces  nouveaux  materiaux  devront  satisfaire  a  des  normes  de  vulnErabilitE  toujours  plus 
strictes.  Qu i  plus  est,  pour  assurer  un  entrainement  soutenu  de  nos  troupes,  ces  nouvelles 
molecules  devront  etre  non  toxiques  et  etre  sans  danger  pour  Penvironnement. 

De  tels  defis  sont  difficiles  a  relever,  car  ces  exigences  semblent  parfois  contradictoires.  Par 
exemple,  on  souhaite  une  plus  grande  performance  tout  en  desirant  une  sensibilite  moindre. 
Pour  repondre  a  ce  defi,  on  devra  creer  des  outils  qui  allegeront  le  travail  du  concepteur  et  lui 
permettront  de  choisir  parmi  plusieurs  candidats.  En  plus  d’etre  couteuse  et  fastidieuse,  la 
conception  de  nouveaux  materiaux  energetiques  produit  une  grande  quantite  de  dechets  qui 
devront  etre  EliminEs  en  respectant  des  normes  environnementales  de  plus  en  plus  strictes  et 
onEreuses.  De  surcroit,  ces  efforts  pourraient  s’avErer  inutiles  si  l’on  constatait  que  les 
produits  manufactures  ne  possEdaient  pas  les  caractEristiques  souhaitEes  ou  prEvues. 

II  serait  done  utile  de  pouvoir  Evaluer  a  priori  certaines  propriEtEs  des  matEriaux  EnergEtiques, 
notamment  les  propriEtEs  mEcaniques  des  mElanges  liant-plastifiant  et,  en  particulier,  la 
tempErature  de  transition  vitreuse  :  Tg.  Cette  demiere  permet  de  prEdire  la  gamme  de 
tempErature  dans  laquelle  un  matEriau  EnergEtique  —  un  propergol  ou  un  explosif  plastique 
—  conservera  des  propriEtEs  mecaniques  acceptables.  Predire  la  tempErature  de  la  transition 
vitreuse  permettrait  une  sElection  prEalable  des  produits  les  plus  prometteurs  et  Eviterait  les 
pertes  de  temps  et  d’argent  dEcoulant  de  travaux  sur  des  substances  inadEquates. 

Ce  travail  dEcrit  justement  les  efforts  effectuEs  au  RDDC  Valcartier  en  utilisant  la  chimie 
computationnelle  afin  d’obtenir  a  priori,  les  tempEratures  de  transition  vitreuses  pour  deux 
Hants  extremement  «  populaires  »  pour  la  formulation  de  compositions  EnergEtiques  soit  le 
polybutadiene  a  terminaison  hydroxyle  et  le  polyazoture  de  glycidyle. 

Ces  travaux  innovateurs  ouvrent  la  porte  a  la  simulation  molEculaire  —  un  nouvel  outil  pour 
les  sciences  de  la  dEfense.  Nous  croyons  que  cet  outil  sera  des  plus  utiles  et  pourra  satisfaire 
les  besoins  de  TarmEe  de  demain  dans  des  domaines  aussi  variEs  que  la  puissance  de  feu,  la 
vulnerabilitE  et  la  sEcurite  des  munitions. 


Abou-Rachid,  H.,  Lafleur-Lambert,  X.  Lussier,  L.-S.,  Mounir  Jaidann,  Gagnon,  N.  et  Brisson,  J. 
2008.  Effet  des  plastifiants  sur  la  tempErature  de  transition  vitreuse,  rg,  du  butadiEne  k  terminaison 
hydroxy  et  du  polyazoture  de  glycidyle  —  Simulations  et  expEriences.  TR  2007-615  RDDC 
Valcartier. 
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1.  Introduction 


Les  prochaines  generations  d’explosifs  a  liant  plastique  devront  presenter  une  plus  grande 
insensibilite  et  une  plus  grande  densite  energetique  tout  en  conservant  une  bonne  integrite 
mecanique. 

En  recourant  a  la  moderation  atomistique  par  la  dynamique  moleculaire  de  systemes 
amorphes  on  peut  predire  les  proprietes  importantes  de  ces  materiaux  ,  et  ce,  avant  leur 
synthese.  On  peut  accelerer  le  travail  de  conception  de  tels  explosifs  en  predisant,  entre  autres, 
Peflfet  d’un  plastifiant  sur  le  polymere  de  base  de  la  composition.  Puisque  les  proprietes  qui 
nous  interessent  dependent  surtout  des  forces  de  dispersion  intermoleculaire,  on  peut  les 
modeliser  &  Paide  de  m^thodes  fondees  sur  la  mecanique  moleculaire. 

L’un  des  objectifs  de  la  mise  au  point  de  nouveaux  explosifs  a  liant  plastique  (PBX)  est  le 
remplacement  du  polymere  inerte  actuellement  utilise,  le  polybutadiene,  par  un  polymere 
energetique  comportant  des  groupements  energetiques,  comme  le  nitrate,  -0N02  ou  Pazoture, 
-N3.  Toutefois,  les  proprietes  elastomeriques  de  ces  polymeres  energetiques  sont  inferieures  a 
celles  du  polybutadiene,  ce  qui  accroit  leur  sensibilite.  On  peut  remedier  a  ce  probleme  en 
ajoutant  une  certaine  proportion  de  plastifiant  au  melange.  On  vise  une  baisse  de  la  valeur  de 
7g,  car  un  materiau  plus  caoutchouteux  absorbe  mieux  les  chocs,  ce  qui  diminue  les  risques  de 
detonation  accidentelle.  On  choisira  un  plastifiant  en  tenant  compte  simultaneities  de  sa 
miscibilite  avec  le  polymere  et  de  sa  capacite  a  baisser  la  valeur  de  la  temperature  de 
transition  vitreuse  d’un  materiau.  La  miscibilite  du  polymere  liquide  et  du  plastifiant  est 
essentielle  a  P integrite  du  materiau  energetique.  En  effet,  puisque  le  melange  de  polymere  et 
de  plastifiant  est  liquide  avant  son  «  durcissement  »  ( curing ),  il  enrobe  completement 
Poxydant  solide  et  les  cristaux  energetiques,  ce  qui  empeche  la  formation  de  vides  dans  le 
materiau. 

Nous  avons  modelise  deux  polymeres,  le  polybutadiene  et  le  polyazoture  de  glycidyle  (connu 
plus  sous  Pacronyme  en  anglais  GAP  au  lieu  de  PAG),  sous  forme  pure  et  dans  un  melange 
50:50  en  masse  avec  un  plastifiant  qui  leur  est  respectivement  compatible.  Dans  notre 
simulation,  le  polybutadiene  est  melange  a  adipate  de  bis(ethyl-2  hexyl)  (DOA),  et 
polyazoture  de  glycidyle  a  du  Gpl,  une  unite  oligomere  du  GAP.  Leur  structure  chimique  est 
montree  &  la  figure  1.  Nous  comparons  les  Tg  predites  avec  les  des  mesures  experimentales 
obtenues  par  analyse  calorimetrique  differentielle  a  compensation  de  puissance  (DSC). 
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Polyazoture  de  glyeidyle  (GAP  et  Gpl) 


Adipate  de  bis(ethyl-2  hexyl) 


Figure  1.  Les  composes  modelises 
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2.  Methodologie 


2.1  Champ  de  forces  et  energie  minimum 

La  modelisation  moleculaire  atomistique  est  une  excellente  methode  de  prediction  de 
differentes  proprietes  des  materiaux,  notamment  la  temperature  de  transition  vitreuse  des 
polymeres.  Cette  modelisation  tient  compte  de  la  position  des  noyaux  atomiques  et  emploie 
un  champ  de  force,  constitue  d’une  serie  d’equations  de  mecanique  classique  pour  estimer  les 
forces  agissant  sur  les  noyaux.  Un  type  est  assigne  a  chaque  atome  auquel  sont  associes  les 
parametres  necessaires  au  calcul  de  1’ energie  par  les  equations  du  champ  de  force.  Celui-ci 
tient  partiellement  compte  de  la  presence  des  electrons,  d’abord  par  l’attribution  d’une  charge 
a  chaque  atome  en  fonction  du  groupement  chimique  auquel  il  appartient,  mais  aussi  en 
fonction  du  type  d’atome.  Souvent,  plus  de  vingt  types  d’atomes  de  carbones  sont  ainsi 
presents,  ce  qui  permet  de  distinguer  les  atomes  de  carbone  sp3,  sp2  ou  sp,  appartenant  a  des 
groupements  aliphatiques  ou  lies  a  des  atomes  d’oxygene,  de  soufre,  d’azote,  etc.  Par 
exemple,  un  atome  de  carbone  specifique  sera  attribue  a  un  groupement  amide,  un  autre  pour 
les  CH2  d’un  cycle  aromatique,  etc. 

Nous  avons  utilise  le  champ  de  force  COMPASS,  offert  dans  le  logiciel  Materials  Studio  4.0 
de  Accelrys  Inc.,  puisqu’il  a  ete  specifiquement  con^u  pour  la  modelisation  moleculaire  de 
materiaux  a  l’etat  solide  [1].  D’ailleurs,  d’autres  auteurs  ont  utilise  ce  meme  champ  de  force 
pour  calculer  des  proprietes,  dont  la  temperature  de  transition  vitreuse,  a  pression,  volume  et 
temperature  constants  de  polyphosphazenes  [2],  de  la  polyisoprene  [3]  et  du  polysilanes  [4]. 

Ce  champ  de  force  dit  de  deuxieme  generation  a  ete  parametrise  a  partir  de  donnees 
experimentales  et  de  resultats  de  calculs  ab  initio  sur  une  grande  variete  de  molecules 
contenant  des  ions  H,  C,  N,  O,  S,  P,  halogenes  et  certains  cations  divalents.  11  a  ete  elabore  a 
partir  du  champ  de  force  PCFF,  lui-meme  issu  du  champ  CFF91.  Chacun  de  ces  champs  de 
force  contient  les  memes  equations  de  bases,  reproduites  a  la  figure  2,  mais  ils  different  par 
leur  parametrisation.  Les  equations  des  divers  termes  du  champ  de  force  sont  donnees  a 
P annexe  2. 

Une  fois  une  molecule  creee  a  l’ecran,  rien  ne  garantit  evidemment  que  son  energie  est 
minimale.  On  doit  done  utiliser  des  methodes  de  minimisation  classiques  pour  d^couvrir  le 
veritable  minimum  de  la  fonction  d’energie.  Puisque  Ton  ne  dispose  d’aucune  methode 
permettant  de  trouver  le  minimum  absolu  de  fonctions  aussi  complexes,  on  se  toume  vers  les 
methodes  mathematiques  iteratives  conduisant  au  minimum  local.  Afin  de  maximiser 
l’efficacite  des  calculs,  on  determine  l’energie  minimale  d’une  molecule  par  une  serie  de 
methodes  sophistiquees.  L’ option  Smart  Minimiser  du  logiciel  Materials  Studio  effectue 
automatiquement  le  choix  de  la  meilleure  sequence  de  methodes. 

Puisque  la  molecule  modelisee  est  normalement  isolee  dans  le  vide  a  0  K,  l’energie  minimale 
calculee  n’est  pas  representative  de  la  realite,  done  les  resultats  de  l’approche  utilisee 
demeurent  qualitatifs. 
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2.2  Dynamique  moleculaire  atomistique 


Les  methodes  de  dynamique  moleculaire  ont  ete  proposees  pour  simuler  les  mouvements  de  la 
molecule  a  une  temperature  donn£e.  Dans  ces  methodes,  on  attribue  arbitrairement  une  vitesse 
initiale  a  chaque  atome.  On  utilise  ensuite  les  equations  classiques  de  Newton  pour  simuler  le 
mouvement  subsequent  des  atomes  sur  la  surface  d’energie  potentielle,  definie  par  le  champ 
de  force. 

Le  calcul  essentiel  de  ces  simulations  est  la  resolution  de  l’equation  de  mouvement  de 
Newton  : 


Fl(t)  =  m,al(t) 


(1) 


Ou  F )  est  la  force  exercee  sur  1’atome  /,  m,  sa  masse  et  al  l’acceleration  qu’il  subit.  La  force 
peut  etre  obtenue  par  la  derivee  de  1’energie  potentielle,  telle  que  definie  par  le  champ  de 
force,  par  rapport  aux  coordonnees  rt  de  1’ atome. 


dV 


(2) 


D’un  point  de  vue  mathematique,  la  principale  difficult^  de  ces  methodes  est  I’impossibilite 
de  calculer  exactement  la  derivee,  etant  donne  la  complexity  de  1 ’equation  de  I’energie.  On 
resout  habituellement  ce  probleme,  en  utilisant  l’equation  de  Verlet  qui  calcule  la  position  et 
la  vitesse  des  atomes  a  un  temps  /  +  A/,  a  partir  de  leur  valeur  au  temps  /.  On  doit  toutefois 
supposer  que  la  vitesse  reste  constante  durant  l’intervalle  A/,  ce  qui  impose  un  intervalle  tres 
court,  de  l’ordre  de  la  femtoseconde. 


II  est  possible  de  fixer  certaines  variables  lors  des  analyses  et,  ainsi,  creer  des  ensembles 
thermodynamiques  distincts.  Le  plus  courant  est  F ensemble  NVE,  pour  lequel  le  nombre 
d’atomes  N,  le  volume  V  et  fenergie  E  sont  maintenus  constants.  C’est  la  methode  la  plus 
rapide  mathematiquement,  mais  la  moins  utile  pour  la  simulation  des  systemes  reels. 
L’ensemble  NVT,  dans  lequel  le  nombre  d’atomes,  le  volume  et  la  temperature  sont  tenus 
constants,  —  aussi  appele  ensemble  canonique  —  exige  une  equation  supplemental  servant 
a  maintenir  constante  la  temperature  (qui  depend  de  la  vitesse  des  atomes). 

L’ensemble  NPT  est  parfois  utilise  pour  modeliser  un  systeme  condense,  ou  le  nombre 
d'atomes,  la  pression  et  la  temperature  sont  tenus  constants.  Le  volume  du  systeme  varie,  ce 
qui  permet  de  changer  la  pression  sur  le  systeme.  Le  controle  de  la  pression  necessite  des 
approches  mathematiques  speciales,  dont  les  deux  principals  sont  dues  a  Berendsen  [5] 
Andersen  [6]. 
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2.3  Simulation  d’un  systeme  amorphe  pur  ou  d’un  melange 


On  simule  la  phase  amorphe  a  l’aide  de  la  methode  de  Theodorou  [7]  qui,  le  premier,  a 
propose  qu’on  puisse  modeliser  une  telle  phase  a  partir  d’une  maille  —  le  plus  souvent 
cubique  —  contenant  des  chaines  disposees  aleatoirement  et  a  laquelle  on  impose  des 
conditions  aux  limites  periodiques.  A  chaque  atome  sortant  de  la  maille  correspond  un  atome- 
image  qui  y  penetrera  par  le  cote  oppose.  On  reproduit  ainsi  l’idee,  imitee  de  la  phase 
cristalline,  qu’une  maille  elementaire  (la  maille  cristallographique)  puisse,  en  etant  repetee  a 
l’infini,  reproduire  un  materiau.  Toutefois,  dans  le  cas  de  la  phase  amorphe,  la  maille  doit  etre 
suffisamment  grande  pour  reproduire  le  phenomene  souhaite,  or  la  dimension  de  cette  «  taille 
suffisante  »  fait  encore  l’objet  de  debats  relativement  aux  temps  des  calculs  de  simulation 
necessaire  pour  atteindre  le  seuil  de  convergence  de  l’energie  minimale.  Que  1  que  soit  le  cas, 
on  introduit  une  symetrie  artificielle,  que  Ton  compense,  notamment,  en  s’assurant  que  la 
demi-hauteur  de  l’arete  de  la  maille  est  superieure  a  la  grandeur  du  phenomene  que  Ton  veut 
reproduire.  A  cause  des  grandes  economies  le  temps  de  calcul  qu’elle  permet,  elle  est  la 
methode  prefer^e  de  la  simulation  des  phases  amorphes.  11  importe  cependant  d’assigner  la 
maille  modelisee  une  dimension  telle  que  les  effets  de  repetition  se  manifestent  dans  un  autre 
ordre  de  grandeur  que  celui  du  phenomene  que  Ton  desire  calculer. 

La  construction  de  la  maille  de  simulation  de  la  phase  amorphe  se  fait  done  en  deux  etapes  :  la 
creation  de  la  maille  et  la  relaxation  des  chaines.  Afin  de  representer  adequatement  la 
distribution  des  conformations  des  chaines,  Theodorou  et  coll.  [7]  ont  propose  d’ajouter  les 
monomeres  une  liaison  chimique  a  la  fois,  sur  une  ou  plusieurs  chaines  en  croissance  dans  la 
maille,  en  fixant  la  conformation  pour  quelle  respecte  simultaneities  une  probability  fonction 
de  la  thyorie  des  ytats  rotationnels  de  Flory  (ou  R1S,  reprysentative  de  l’energie  interne  de  la 
chaine),  des  contacts  entre  la  chaine  en  croissance,  des  atomes-images  et,  finalement,  de 
l’energie  inter-chaine.  Cependant,  cette  croissance  provoque  parfois  des  recouvrements 
importants  d’atomes  lors  de  la  creation  d’un  systeme  a  partir  d’une  longue  chaine  polymere. 


La  methode  de  Theodorou  construit  les  chaines  dans  une  maille,  aleatoirement  (Monte  Carlo), 
segment  par  segment.  Meirovitch  [8]  a  propose  que  Ton  tienne  compte,  non  seulement  de 
l’energie  d’un  atome  ou  groupe  d’atomes  ajoute  lors  de  la  croissance  de  la  chaine,  mais  aussi 
de  l’energie  de  quelques  groupes  additionnels,  ce  qui  evite  les  situations  menant  a  des  culs-de- 
sac  energytiques,  e’est-a-dire  a  des  solutions  d’energie  anormalement  elevee.  Meirovitch  a 
montre  que  sa  methode  permettait  d’estimer  la  fonction  de  partition  d’un  systeme  amorphe. 

Le  module  Amorphous  Cell  du  logiciel  Materials  Studio  incorpore  des  elements  de  ces  deux 
methodes. 

Ces  methodes  se  pretent  aussi  a  la  fabrication  de  modeles  amorphes  de  melanges.  Dans  tous 
les  cas,  on  obtient  un  systeme  plus  realiste  en  recourant  a  la  dynamique  moleculaire  pour 
relaxer  les  conformations  de  hautes  energies.  On  peut  egalement  reajuster  la  densite  lors  de 
cette  etape. 
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Lors  des  experiences  de  dynamique  moleculaire,  P  ensemble  adiabatique  NVE  n’est  utile  que 
dans  de  rares  cas,  par  exemple  si  Ton  souhaite  explorer  la  surface  d’energie  constante  de 
Pespace  conformationnel  du  systeme.  Cet  ensemble  ne  permet  cependant  pas  d’equilibrer  un 
systeme,  car  il  n’y  a  pas  d’echange  d’energie  avec  Texterieur  et,  ainsi,  on  ne  peut  pas  obtenir 
une  density  specifique  correspondant  a  une  temperature  donnee.  L’ensemble  NVT  permet 
d’obtenir  l’equilibre  energetique  d’un  systeme,  mais  pas  sa  relaxation  complete,  car  le  volume 
est  garde  constant  (done  sa  densite).  11  est  utile  lorsqu’on  souhaite  maintenir  la  densite  a  une 
valeur  donnee.  L’ensemble  NPT  est  le  plus  frequemment  utilise  pour  les  simulations  :  il 
permet  effectivement  a  un  systeme  ayant  une  temperature  donnee  d’atteindre  la  densite 
d’equilibre,  tout  en  conservant  T  fixe  durant  toute  la  simulation.  Toutefois,  on  se  souviendra 
que  la  densite  varie  brusquement  lors  de  transition,  y  compris  la  transition  vitreuse  dont  nous 
voulons  estimer  la  temperature  pour  differents  melanges,  avant  la  synthese  en  laboratoire  des 
meilleurs  candidats  pour  certaines  applications. 


2.4  Determination  de  Tg  par  la  dynamique  moleculaire 


Les  transitions  thermodynamiques  imposent  des  modifications  a  diverses  proprietes, 
notamment  la  temperature  de  fusion  —  une  transition  du  premier  ordre.  La  transition  vitreuse 
est  assimiiee  a  une  transition  du  second  ordre.  On  a  propose  de  nombreux  modeles  pour 
expliquer  la  temperature  de  transition  vitreuse.  Dans  certains,  la  transition  vitreuse  est 
presentee  comme  un  phenomene  thermodynamique  [9];  dans  d’autres  comme  un  processus 
cinetique  [10].  Le  modele  du  volume  libre  [11],  combine  adroitement  les  composantes 
thermodynamiques  et  cinetiques  de  cette  transition,  mais  une  discontinuite  des  proprietes 
d’etat  en  fonction  de  la  temperature  pourrait  apparaitre.  L’effet  cinetique  est  clairement  etabli 
et  superpose  a  une  telle  discontinuite  [12]. 

Rigby  et  Roe  ont  demontre  que  les  modeles  du  polyethylene  obtenus  par  modeiisation 
moleculaire  presentaient  une  telle  discontinuite  pour  le  volume  specifique  et  les  coefficients 
de  diffusion  [13],  ainsi  que  pour  la  distribution  du  volume  libre  [14].  Par  la  suite,  plusieurs 
chercheurs  ont  demontre  que  ces  discontinuites  apparaissaient  egalement  chez  d’autres 
polymeres  [15-17].  Boyd  a  publie  une  synthese  de  la  litterature  sur  ces  questions  [18].  Ces 
conclusions  ont  ete  contestees,  puisque  la  vitesse  de  refroidissement,  dans  les  experiences  de 
modeiisation,  est  nettement  inferieure  au  temps  de  relaxation,  qui  est  de  l’ordre  de  10'6  s. 

Dans  ce  contexte,  il  est  clair  que  Petat  d’equilibre  ne  peut  etre  atteint  et  que  les  proprietes 
calcuiees  ne  peuvent  representer  celles  qui  existant  a  l’equilibre  [19].  D’une  part,  en  depit  ces 
critiques,  au  fur  et  a  mesure  des  travaux,  il  apparait  que  la  modeiisation  moleculaire  simule 
suffisamment  bien  les  mouvements  responsables  ou  resultants  de  P  augmentation  du  volume 
libre  a  la  Tg,  pour  croire  que  les  resultats  de  ces  calculs  sont  une  bonne  estimation  de  la  realite. 
D’autre  part,  tous  s’accordent  sur  le  fait  qu’il  existe  inevitablement  un  ecart  variable  entre  la 
temperature  de  transition  vitreuse  reelle  et  le  calcul  de  celle-ci.  D’autres  auteurs  ont  applique 
cette  technique  de  modeiisation  a  des  melanges  [20]  et,  recemment,  a  la  prediction  de  la  T% 
d’un  polymere  en  presence  de  plastifiant  [21]. 

Plusieurs  conditions  peuvent  modifier  la  temperature  de  transition  vitreuse.  Toutefois,  les 
systemes  liquides  et  solides  ne  sont  affectes  que  par  des  pressions  de  l’ordre  du  gigapascal 
GPa,  or  la  pression  atmospherique  normale  se  situe  autour  de  0,0001  GPa.  C’est  pourquoi, 
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nous  avons  attribue  une  pression  exteme  de  0  GPa  a  T  ensemble  NPT,  qui  pourrait  etre 
consideree  comme  une  approximation  acceptable  des  conditions  normales  [1]  puisque  les 
faibles  variations  de  pression  qui  y  prevalent  n’affecteront  pas  les  resultats.  Le  temps  de 
simulation  influence  la  valeur  de  la  Tt,  mais  on  a  montre  que  Teffet  de  cette  variable  etait 
minime  [2].  La  valeur  de  Tg  obtenue  depend  surtout  de  Tequilibre  du  modele  de  depart  :  un 
modele  mal  equilibre,  ou  une  maille  trop  petite  donneront  des  valeurs  qui  divergent  davantage 
des  valeurs  exp^rimentales.  Fried  et  Ren  ont  demontre  que  les  valeurs  d’expansion  thermiques 
etaient  tres  sensibles  a  la  taille  de  la  maille,  et  que  la  Tg  des  polysilanes  etait  surestimee 
d’environ  15  K  lorsque  la  maille  etait  trop  petite.  11s  ont  propose  que  des  mailles  d’environ  35 
a  40  devraient  permettre  d’obtenir  des  valeurs  de  Tg  fiables.  Lors  de  la  construction  des 
systemes,  on  doit  done  s’assurer  que  la  taille  des  modeles  soit  egale  a  cette  valeur  ou  la 
depasse. 

Pour  effectuer  les  calculs  de  transition  vitreuse  par  modelisation  moleculaire,  on  doit  d’abord 
effectuer  un  second  equilibrage  de  100  a  150  ps,  a  100  K  au-dessus  de  la  valeur 
approximative  de  Tg,  avec  les  memes  options  que  precedemment.  La  temperature  sera 
diminuee  par  pas  de  10  K,  toujours  en  NPT,  et  on  reequilibrera  chaque  temperature  durant 
150  ps.  La  densite  du  systeme  sera  alors  obtenue  par  la  moyenne  sur  les  50  demieres 
picosecondes  de  la  simulation.  On  pourra  alors  tracer  une  courbe  de  volume  sp^cifique  en 
fonction  de  la  temperature. 


2.5  Procedures  de  calcul  de  modelisation 


Tous  les  calculs  ont  ete  effectues  avec  le  logiciel  Materials  Studio  de  Accelrys  Inc.,  installe 
sur  un  ordinateur  SGI  Altix  3000  du  Centre  de  calcul  bio-informatique  de  TUniversite  Laval. 
Le  champ  de  force  COMPASS  a  ete  utilise  dans  tous  les  calculs. 


2.5.1.  Creation  des  chames  polymeres 


On  doit  definir  l’unite  de  repetition  (Build  .  Build  Polymers  :  Repeat  Unit),  si  elle  n’est  pas 
dejei  precisee  dans  le  logiciel.  (Voir  le  protocole  des  calculs  a  Tannexe  2.)  On  doit  definir  les 
des  atomes  de  tete  ( head  atom)  et  de  queue  ( tail  atom)  de  V unite  de  repetition,  qui 
repr^sentent  les  points  de  liaisons  des  monomeres  pour  assembler  le  futur  polymere.  Ces 
atomes  terminaux  doivent  etre  des  hydrogenes.  Lors  de  la  polymerisation,  ils  sont  retires  pour 
faire  place  h  une  liaison  avec  le  monomere  suivant. 

Le  polymere  est  ensuite  construit  par  polymerisation  (Build  Build  Polymers  :  Homopolymer) 
en  specifiant  l’unite  de  repetition  (Unit :  Current  Project)  et  le  nombre  d’unites  de  repetition. 
Si  on  desire  placer  un  atome  autre  que  1’hydrogene,  on  doit  specifier  la  tacticite  et  les 
groupements  initiateur  et  terminal  de  la  chaine.  On  peut  aussi  fixer  Tangle  de  torsion  de  la 
liaison  creee  entre  les  monomeres.  11  n’est  pas  recommander  de  specifier  un  angle  de  torsion 
autre  que  180°,  dans  le  cas  des  polymeres  possedant  des  groupes  lateraux,  on  risque  de 
produire  une  molecule  dont  Tenergie  est  anormalement  eievee.  Chaque  structure  sera 
minimisee  a  son  tour,  ce  qui  la  placera  dans  une  conformation  stable. 
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Une  fois  le  copolymere  ou  un  homopolymere  construit,  avec  une  tacticite  donn£e,  on  passe  a 
l’6tape  d’assemblage  des  copolymeres  (Build :  Build  Polymers  :  Random  Copolymer),  en 
specifiant  les  unites  de  repetition  (qui  peuvent  etre  les  unites  S  et  R  dans  le  cas  de  variations 
de  tacticite)  et  en  leur  assignant  les  probability  desirees. 

Une  fois  le  polymere  cree,  on  le  soumet  a  une  minimisation  ( Discover  :  Minimizer)  par  la 
methode  Smart  Minimizer  et  une  convergence  ultrafine.  Cette  option  permet  d’effectuer  la 
minimalisation,  en  employant  des  methodes  mathematiques  de  plus  en  plus  precises,  et  ce,  de 
fagon  automatis^e. 

Cette  procedure  de  creation  et  de  minimisation  du  systeme  est  r&teree  a  trois  reprises,  afin  de 
creer  au  moins  trois  chaines  optimisees.  On  ne  conserve  que  l’essai  ayant  donn6  Penergie  la 
plus  basse.  Ceci  permet  fr  garder  seulement  les  conformations  du  systeme  ayant  converge  vers 
un  minimum  local  d’energie  et  d’eliminer  les  autres.  Une  telle  approche  permettrait  d’avoir 
des  conformations  plus  representatives  de  la  physique  des  systemes  reels. 

Des  chaines  de  polybutadiene  a  radical  hydroxy  ont  ete  construites  par  Paddition  aieatoire  de 
60  %  d’unites  trans-1,4,  20%  d’unites  cis-1,4  et  20%  d’unites  vinyl-1, 2  .  Le  nombre  d’unites 
de  monomeres,  ainsi  que  la  masse  des  molecules  construites  sont  indiques  au  tableau  1.  Une 
fois  la  chaine  creee,  on  la  soumet  a  une  minimisation  (Discover  :  Minimizer)  avec  la  methode 
Smart  Minimizer  a  un  niveau  de  convergence  ultra-fine.  Trois  chaines  ont  ainsi  ete  creees  et 
minimisees. 


2.5.2  Procedure  de  creation  de  modeles  de  la  phase  amorphe 


Pour  creer  un  modele  de  phase  amorphe,  on  doit  d’abord  choisir  une  chaine  polymere  dans 
une  conformation  de  basse  energie  qui  servira  a  la  creation  de  la  cellule  amorphe  (Amorphous 
Cell :  Construction).  On  specific  le  nombre  de  molecules  pour  obtenir  des  cellules  amorphes 
contenant  au  moins  1000  atomes  environ,  puisque  le  champ  de  force  COMPASS  est  etalonne 
pour  des  systemes  de  cette  taille.  On  doit  egalement  creer  au  moins  trois  modeles  differents 
(ou  trois  configurations)  des  phases  amorphes  —  idealement  9  ou  10  —  pour  s’assurer  d’avoir 
un  systeme  dont  l’energie  serait  anormalement  eievee.  (Une  consequence  de  Palgorithme  de 
construction  aieatoire  de  cellules  amorphes).  On  minimise  Penergie  de  chacun  des  modeles 
construits,  et  on  choisit  ensuite  le  systeme  de  plus  basse  energie  pour  la  relaxation  par 
dynamique  moieculaire.  Normalement,  si  elle  est  connue,  on  foumit  au  programme  pour  la 
construction,  la  masse  volumique  determinee  experimentalement.  S’il  s’agit  d’un  melange,  on 
utilisera  la  moyenne  pond^ree  des  masses  volumiques  de  chaque  substance.  Si  Pon  ne  connait 
pas  la  masse  volumique  de  la  substance,  on  peut  utiliser  la  valeur  par  defaut  de  1  g/cnr  .  Apres 
a  la  construction,  on  procedera  a  une  simulation  par  dynamique  moieculaire,  dans  le  sous- 
ensemble  NPT,  pour  laisser  le  volume  converger  vers  sa  valeur  d’equilibre.  Habituellement, 
cette  valeur  reproduit  assez  bien  les  donnees  experimentales  lorsque  Pon  utilise  un  champ  de 
force  appro prie  pour  la  matiere  condensee,  ce  qui  est  le  cas  du  champ  de  force  COMPASS  [1], 

Une  fois  les  modeles  crees,  on  doit  s’assurer  que  les  molecules  sont  bien  relaxees.  On  suivra 
done  Involution  Penergie  du  systeme  et  de  son  volume,  jusqu’a  ce  que  la  masse  volumique  et 
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l’energie  potentielle  restent,  en  moyenne,  constantes  durant  le  calcul  de  la  dynamique 
moleculaire  dans  l’ensemble  NPT. 

La  temperature  est  maintenue  constante  par  la  methode  d’ Anderson,  et  la  pression  par  la 
methode  de  Berendsen.  Puisque  nous  simulons  un  systeme  infini,  avec  conditions  des  limites 
periodiques,  il  est  necessaire  d’etablir  une  distance  pour  les  calculs.  Nous  utilisons 
generalement  les  options  par  defaut  du  logiciel,  done  une  distance  limite  de  9,5  A  pour  les 
interactions  non  liees  de  type  van  der  Waals,  apres  laquelle  nous  appliquons  une  fonction 
spline  degressive  d’une  largeur  de  0,5  A  sur  le  1,0  A  suivant.  Le  calcul  des  interactions 
electrostatiques,  qui  se  font  sentir  a  plus  longue  distance,  est  fait  par  sommation  d’Ewald 
(dans  l’espace  reciproque,  par  transformee  de  Fourier),  avec  une  precision  de  10  2  et  une 
largeur  de  mise  a  jour  de  1,00.  Cette  methode  permet  d’evaluer  avec  une  grande  precision  la 
valeur  qu’aurait  un  systeme  non  fini.  Toutefois,  l’utilisation  de  la  methode  d’Ewald  fait 
exploser  le  temps  de  calcul  aussi,  pour  de  plus  gros  systemes,  il  est  parfois  necessaire  de 
proceder  une  approximation  type  employee  pour  les  energies  de  Van  der  Waals.  L’integration 
se  fait  selon  l’algorithme  de  Verlet,  avec  un  pas  d’une  femtoseconde. 

Le  protocole  de  creation  de  la  maille  amorphe  est  presente  en  annexe.  On  peut  ajuster  les 
parametres  suivants  : 

•  Maximum  lookahead  bonds  (nombre  maximum  de  liaisons  considerees  lors  de 
l’etablissement  de  la  conformation  d’un  monomere  ajoute) :  La  valeur  par  defaut 
est  1,  ce  qui  correspond  a  la  methode  classique  de  Theodorou  et  Suter.  On 
l’augmenter  jusqu’a  6,  si  Ton  suit  la  methode  de  Meirovitch. 

•  Maximum  lookahead  configurations  (nombre  maximum  de  conformations 
considerees  lors  de  l’ajout  d’un  monomere):  Par  defaut,  cette  valeur  est  de  48,  mais 
elle  peut  varier  de  1  a  128.  Elle  correspond  au  nombre  de  conformations  que  le 
logiciel  peut  essayer,  chaque  fois  qu’il  additionne  un  monomere,  avant  ne  plus 
pouvoir  en  ajouter. 

•  Number  of  bonds  to  be  added  per  step  (nombre  de  liaisons  a  ajouter  a  chaque  etape) : 

Il  permet  de  specifier  le  nombre  d’etapes  a  effectuees  lors  de  la  construction.  On 
ajoute,  par  defaut,  une  liaison  chimique  a  la  fois. 

•  Number  of  substates  per  state  (nombre  de  sous-etats  conformationnels)  :  La  valeur 
par  defaut  est  1,  mais  on  peut  1’accroitre  jusqu’a  15. 

•  Substate  width  (largeur  de  sous-etat) :  Nombre  de  degres  de  chaque  sous-£tat.  Par 
defaut,  20,  mais  cette  valeur  peut  varier  de  1  a  1  80. 

•  Check  for  ring  catenations  (verification  de  la  catenation  en  cycle)  Ce  parametre 
demande  au  programme  de  verifier,  lors  de  l’ajout  d’un  monomere,  si  celui-ci  est 
passe  au  milieu  d’un  cycle,  ce  qui  conduit  a  des  etats  d’energie  Aleves,  improbables 
dans  la  nature,  mais  qu’il  sera  difficile  d’^liminer,  par  dynamique  moleculaire.  dans  la 
structure  finale. 

Nous  specifions  le  nombre  de  molecules  (ou  chaines  de  polymeres,  selon  le  cas).  11  est 
conseille  de  construire  les  phases  amorphes  avec  une  centaine  d’atomes  environ,  puisque  le 
champ  de  force  COMPASS  a  ete  etalonne  pour  des  systemes  de  cette  taille.  De  plus,  des 
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systemes  trop  grands  exigent  des  temps  de  calcul  enormes  et  la  memoire  necessaire  peut 
d^passer  les  capacites  de  Fordinateur  ou  du  logiciel. 

Sur  une  base  de  construction  arbitraire  tenant  compte  de  F  existence  de  toutes  les 
conformations  dans  une  maille,  le  polymere  polybutadiene  est  simule  sous  la  forme  d’un 
melange  de  trois  unites  de  repetition  differentes  :  60%  dans  la  configuration  trans-1,4  qui  est 
generalement  la  configuration  la  plus  stable  et  probable,  20%  dans  la  configuration  cis-1 ,4 
moins  stable  en  general  a  cause  des  effets  steriques  et  20%  dans  la  configuration  vinyl-1 ,2  qui 
est  aussi  aussi  probable  que  la  cis.  Le  polybutadiene  est  alors  construit  en  utilisant  la 
commande  :  Build  :  Build  Polymers  :  Random  Copolymer  qui  servent  a  specifier  les  trois 
unites  de  repetition  et  en  leur  assignant  les  probability  0,6  pour  la  configuration  trans,  0,2 
pour  la  cis  et  0,2  pour  la  vinyl.  Une  fois  le  polymere  cree,  on  la  soumet  a  un  calcul  de 
minimisation  avec  la  commande  Discover  :  Mmimizer  en  cochant  les  options  Smart 
Minimizer  et  le  niveau  de  convergence  Ultra-fine.  Cette  procedure  de  creation  et  de 
minimisation  est  reiter^e  trois  fois  afin  de  cr£er  au  moins  trois  chaines  optimisees. 

Egalement,  il  est  aussi  essentiel  de  donner  une  masse  volumique  initiale  au  systeme  qu’on 
desire  creer.  Generalement,  on  lui  attribuera  la  masse  volumique  experimentale  si  elle  est 
connue.  Si  elle  etait  inconnue,  on  peut  utiliser  la  valeur  par  defaut  de  1  g/cm3.  A  la  fin  des 
calculs,  on  utilisera  la  dynamique  moleculaire  pour  donner  une  valeur  plus  realiste  a  cette 
masse  volumique. 

Une  fois  les  mailles  construites,  il  est  conseille  de  choisir  celle  de  plus  basse  energie  pour  le 
reste  des  calculs.  Habituellement  avec  l’ensemble  NPT,  on  effectue  une  relaxation  du 
systeme,  avec  des  temps  100  a  150  ps  —  parfois  plus  longtemps  — ,  jusqu’a  ce  que  les 
propri&es  moyennes  du  systeme  (Energie  potentielle,  Energie  exteme  ou  interne,  rayon  de 
giration,  etc.)  ne  varient  plus  de  fa^on  appreciable. 


2.5.3  Procedure  d’equilibration  du  volume  a  differentes 
temperatures 


Une  fois  la  maille  construite  et  equilibree,  des  calculs  de  la  dynamique  moleculaire  sont 
effectors  en  utilisant  Pensemble  NPT  pour  les  differentes  temperatures.  Cet  ensemble  permet 
la  variation  de  Fenergie  du  systeme  et  de  son  volume.  Une  densite  optimale  du  systeme  peut 
ainsi  etre  atteinte  a  chaque  temperature.  La  pression  exteme  est  fixee  0,0  GPa  avec  l’ensemble 
NPT,  ce  qui  constitue  une  simulation  acceptable  des  conditions  normales  [1].  La  temperature 
et  la  pression  sont  maintenues  constantes  selon  les  methodes  d’Anderson  et  de  Berendsen 
respectivement.  A  chaque  temperature,  la  density  de  chaque  systeme  est  calculee  k  partir  de  la 
moyenne  sur  les  50  demieres  picosecondes  du  calcul  de  dynamique  moleculaire. 
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Tableau  1.  Composition  des  mailles  amorphes  construites 


HTPB 

HTPB-DOA 

GAP 

GAP-Gpl 

Gpl 

Molecules  de  polymere 

2 

1 

1 

1 

0 

Molecules  de  plastifiant 

0 

8 

0 

7 

12 

Nombre  d’unit6s  de  monom^re 
par  chaine  polymere 

53 

53/1 

50 

50/7 

7 

Nombre  total  d’atomes  dans  la 
maille 

1028 

1058 

603 

1205 

1032 

Fraction  massique  de 
plastifiant 

0,000 

0,514 

0,000 

0,495 

1,000 

Fraction  molaire  de  plastifiant 

0,000 

0,875 

0,000 

0,857 

1,000 

Density  experimental  (g/cm3) 

0,899 

— 

1,300 

— 

1,240 

Density  imtiale  (g/cm3) 

0,900 

0,915 

1,300 

1,270 

1,240 

Density  finale  (g/cm3) 

0,893 

0,916 

1,277 

1,288 

1,261 

Parametre  de  maille  final  (A) 

21,83 

21,92 

18,62 

23,40 

22,32 

2.6  Mesures  experimentales 

Nous  avons  mesure  les  temperatures  de  transition  vitreuses  a  l’aide  d’un  calorimetre  Q2000 
de  TA  Instruments.  Chaque  echantillon  contenait  de  4  a  6  mg  de  polymere  pur  ou  d’un 
melange  polymere-plastifiant,  dans  une  cellule  standard  d’analyse  calorimetrique 
differentielle.  Nous  avons  etalonne  Tappareil  en  utilisant  le  point  de  fusion  de  Pindium 
comme  reference.  Chaque  melange  a  ete  elabore  juste  avant  les  mesures,  en  combinant  entre  2 
et  3  mg  de  polymere  avec  une  masse  £gale  du  plastifiant  correspondant.  Nous  avons  effectue 
les  analyses  sous  atmosphere  d’air  purifie  (50  ml/min)  en  elevant  la  temperature  de 
l’echantillon  avec  un  taux  1 0  °C/min  depuis  -1 20  °C  jusqu’ik  1 50  °C.  La  temperature  de 
transition  vitreuse  mesuree  correspond  au  milieu  du  saut  de  la  capacite  calorifique. 

Tableau  2.  Proprietes  des  polymeres  et  plastifiants  utilises 


HTPB 

GAP 

Gpl 

DOA 

Masse  (g/mol) 

2800 

5000 

700 

371 

Polydispersit6 

— 

1,0 

Viscosity  (cps) 

50 

12000 

125 

— 

Fonctionnalit6* 

2,3 

2,7 

0 

— 

*  Nombre  de  groupements  OH  terminaux,  tel  que  specific  par  le  fabricant. 


Le  polybutadiene,  commercialise  sous  le  nom  de  resine  Poly  Bd®  R-45M,  a  ete  fabrique  par 
Sartomer  Inc.,  le  GAP  (polymere)  et  le  Gpl  (plastifiant)  provient  de  3M  Corp.,  et  le  dioctyi 
adipate  est  vendu  par  Canada  Colors  &  Chemicals  Ltd. 
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3.  Res  u  I  tats  et  discussions 


3.1  Determination  experimental  de  la  transition  vitreuse,  et 
methodes  de  prediction  a  partir  de  la  temperature  des 
composes  purs 

Selon  Stacer  et  Husband  [22],  la  valeur  de  la  temperature  de  transition  vitreuse  (Tg)  constitue 
un  critere  de  d^faillance  des  polymeres  du  point  de  leurs  mecaniques  puisque,  sous  cette 
temperature,  le  materiau  perd  rapidement  ses  proprietes  elastomeriques.  Ainsi,  il  est  important 
d’evaluer  la  capacite  de  prediction  de  ce  parametre  et  les  limites  de  precision  de  ces 
predictions,  pour  des  systemes  comportant  des  molecules  energetiques. 

Pour  proceder  a  ces  evaluations,  nous  avons  choisi  deux  materiaux  contenant  des  molecules 
energetiques :  le  melange  polybutadiene-adipate  de  bis(ethyl-2  hexyl)  (HTPB-DOA)  et  le 
melange  Polyazoture  de  glycidyle-monere  (GAP-Gpl),  plus  leger.  Nous  avons  mesure  la 
temperature  de  transition  vitreuse  de  chaque  materiau,  a  l’aide  d’un  calorimetre  analytique  a 
compensation  de  puissance,  dans  des  conditions  bien  controlees,  puisqu’il  est  bien  connu  que 
la  temperature  de  transition  vitreuse  varie  aussi  bien  en  fonction  de  la  masse  mol^culaire  d’un 
polymere  que  la  vitesse  du  balayage  d’analyse  enthalpique  differentielle. 

On  retrouve  dans  la  figure  3  des  thermogrammes  ty piques  des  balayages  d’analyse 
enthalpique  differentielle.  Aucun  balayage  du  DOA  n’est  presente,  puisque  ce  compose  est 
fortement  cristallin  —  sa  phase  amorphe,  dans  laquelle  se  produit  la  transition  vitreuse,  etant 
virtuellement  absente,  aucun  changement  dans  l’absorption  de  chaleur  n’a  ete  observe. 
Comme  on  peut  le  voir  sur  ces  graphiques,  l’ajout  de  DOA  au  HTPB  se  traduit  par 
une  baisse  de  la  temperature  a  laquelle  on  observe  le  saut  de  capacite  calorifique, 
typique  du  phenom&ne  de  transition  vitreuse.  De  plus,  on  note  egalement  que  la 
largeur  de  cette  transition  —  le  domaine  de  temperature  a  Pinterieur  duquel  se  produit 
cette  transition  —  reste  plutot  constante  Ces  deux  observations  indiquent  la  miscibilite 
du  systeme,  du  moins  a  l’^chelle  de  dimensions  de  domaines  de  phases  sdparees  pour 
laquelle  la  calorimetrie  donne  une  seule  transition  de  largeur  constante,  qui  a  etd 
etablie  comme  etant  aux  alentours  de  10-50  nm  [23,  24].  Nous  avons  fait  des 
observations  similaires  sur  le  deuxieme  systeme  :  les  melanges  GAP-Gpl . 

Les  temperatures  de  transition  observees  sont  enumerees  au  tableau  3.  On  remarque  des  ecarts 
de  5  a  10  °C,  par  rapport  aux  donnees  publiees,  qui  pourraient  decouler  des  differences  de 
vitesse  de  balayage;  or,  les  references  n’ indiquent  pas  la  vitesse  utilisee  pour  les  mesures. 
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Figure  3  :  Endothermes  des  composes  purs  et  des  melanges  mesures  a  10°/min  par 
analyse  calorimetrique  a  compensation  de  puissance 
(Heat  Flow  :  transfert  de  chaleur) 


a)  HTPB  (R45M)  pur 


Temperature  (K) 

b)  Melange  HTPB  (R45M)  et  DOA 


Temperature  (  K  ) 
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c)  GAP  pur 


d)  Gpl  pur 
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e)  Melange  GAP  et  GAP  plastifiant 


T emperature  (K) 


On  a  propose  plusieurs  methodes  pour  predire,  la  Tg  d’un  melange  a  partir  des  temperatures  de 
transition  pour  les  substances  pures.  La  plus  simple  est  la  loi  d’additivite  : 


T  =  w.T  .  +  w’T  j 

g  I  g,l  2  g,2 


(3) 


ou  w,  est  la  fraction  massique  du  compose  i  et  Fg<f  sa  temperature  de  transition  vitreuse. 
Plusieurs  autres  equations  ont  ete  proposees,  dont  celle  de  Fox  [25], 


1  =  ,  ™2 
T  T  T 

qui  ne  comporte  aucun  autre  parametre,  et  celle  de  Gordon-Taylor, 

wJ  .  +  kwS 


KiJg,i 


2  g.2 


+  kw2 


(4) 


(5) 


qui  contient  un  terme  correctif,  k ,  tenant  compte  des  changements  du  volume  specifique  et  des 
coefficients  d’expansion  volumique  des  polymeres  a  la  temperature  de  transition  vitreuse.  A 
cause  du  nombre  accru  de  parametre,  Putilisation  de  cette  Equation  est  limitee  aux  predictions 
basees  sur  les  donnees  publiees  [26].  On  a  observe,  chez  les  systemes  presentant  des 
interactions  fortes,  des  ecarts  importants  avec  les  donnees  experimentales,  ce  qui  a  conduit 
Kwei  et  coll.  [27]  a  proposer  un  terme  correctif  a  1 ’equation  de  Gordon-Taylor, 


iv,  +  kw2 


+  ^iv,iv2 


(6) 
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dans  laquelle  q  est  un  terme  ajustable  dependant  de  la  force  des  interactions  entre  les  deux 
composants  du  melange.  Cons6quemment,  seules  la  loi  d’additivite  et  Tequation  de  Fox  sont 
utilisables  lorsque  les  seules  donnees  experimentales  connues  sur  les  composantes  pures  sont 
les  temperatures  de  transition  vitreuses.  Les  valeurs  obtenues  pour  le  seul  systeme  dont  nous 
connaissions  les  temperatures  de  transition  vitreuses,  le  GAP-Gpl,  montrent  un  bon  accord 
avec  la  mesure  experimental.  La  loi  d’additivite  sous-evalue  la  temperature  de  transition 
vitreuse  d’environ  4°,  alors  que  l’equation  de  Fox  la  surevalue  de  2°.  Or,  nous  evaluons  a  5°, 
la  precision  de  notre  mesure  de  Tg ,  en  raison  de  la  largeur  de  la  transition,  de  l’effet  de  la 
vitesse  de  la  mesure  (effet  cinetique),  et  de  la  presence  possible  de  contaminants.  Les  deux 
systemes  choisis  ne  sont  pas  les  meilleurs  candidats  pour  eprouver  la  validite  de  ces 
equations,  puisque  l’ajout  du  plastifiant  modifie  peu  la  temperature  de  transition  vitreuse  du 
melange  (15°  pour  le  systeme  HTPB-DOA  et  14°  pour  le  systeme  GAP-Gpl).  Cet  effet 
modeste  decoule  directement  du  faible  ecart  entre  les  temperatures  de  transition  vitreuse  des 
deux  composants  du  melange  (seulement  22°  entre  le  GAP  et  le  Gpl). 

Cette  methode  rapide  et  efficace  de  prediction  oblige  cependant  de  connaitre  les  temperatures 
de  transition  vitreuse  des  composants  du  melange.  Or,  elles  sont  parfois  inconnues, 
notamment  dans  le  cas  de  petites  molecules  cristallines  comme  le  DOA.  En  outre,  meme  si  les 
temperatures  de  transition  vitreuse  de  deux  composes  etaient  connues,  leur  miscibilite  ne  l’est 
pas  toujours.  De  plus,  une  prediction  pour  un  melange  immiscible,  faite  a  partir  de  ces 
equations  ne  serait  evidemment  pas  valide.  Finalement,  dans  certains  cas,  on  pourra  desirer 
predire  l’effet  de  certaines  molecules  sur  la  transition  vitreuse,  avant  leur  synthese.  Les 
methodes  de  moderation  moieculaire  pourront  alors  etre  utiles,  en  depit  des  longs  temps  de 
calcul. 


Tableau  3.  Comparaisons  entre  les  Tg  obtenues  experimentalement,  les  T%  publiees  dans  les 
ecrits  scientifiques  et  celles  predites  a  partir  des  Tgdes  composes  purs. 


Systeme 

T  exP 
ffl 

K 

T  pub. 

K 

y  add. 

'o 

K 

y  FOX 
'fl 

K 

(±10) 

HTPB 

192 

198  [28] 

— 

— 

HTPB-DOA 

177 

— 

— 

— 

GAP 

227 

237  [4] 

— 

— 

GAP-Gpl 

213 

218  [4] 

209 

215 

Gpl 

205 

— 

— 

— 

3.2  Simulation  des  phases  amorphes  par  la  dynamique 
moieculaire 

La  mecanique  et  la  dynamique  moieculaire  atomistiques  permettent  de  simuler  aussi  bien  une 
phase  amorphe  qu’une  phase  cristalline.  En  outre,  ainsi  que  nous  Fillustrions  dans  un  travail 
anterieur  [29],  on  peut  utiliser  ces  techniques  pour  predire  la  miscibilite  de  molecules 
energetiques. 

L’ utilisation  du  progiciel  Materials  Studio ,  et  plus  precisement  de  son  module  Amorphous 
Phase ,  a  permis  de  cmer  des  modeles  raisonnables  de  la  phase  amorphe.  Un  modele 
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representatif  des  mailles  creees,  le  modele  de  plus  basse  energie,  est  montre  a  la  figure  4  de 
deux  manieres.  Dans  un  premier  temps,  on  peut  voir  les  chaines  polymeres  et  les  petites 
molecules,  telles  que  construites  et  places  dans  la  maille,  mais  on  ne  voit  pas  les  molecules- 
images  produites  lorsque  Ton  applique  les  conditions  limites  periodiques.  Les  molecules 
occupent  peu  de  volume  dans  la  maille  et  la  densite  trouvee  experimentalement  n’est  pas 
reproduite.  On  peut  toutefois  visualiser  la  chaine  ou  la  molecule,  qui  adopte,  lorsque  la  masse 
moleculaire  est  elevee,  une  conformation  sinueuse,  partiellement  repliee  sur  elle-meme,  ce  qui 
est  generalement  le  cas  des  chaines  polymeres.  Dans  la  seconde  representation  du  meme 
modele,  on  voit  les  atomes  des  molecules  de  depart  et  ceux  des  molecules-images  qui 
envahissent  l’interieur  de  la  maille.  L’espace  est  entierement  occupe,  on  ne  voit  aucun  espace 
de  faible  densite.  Si,  dans  le  cas  des  chaines  de  forte  masse  moleculaire,  on  ne  distingue  que 
difficilement  les  chaines  isolees,  avec  de  petites  molecules,  les  deux  representations  sont  tres 
similaires,  les  molecules  debordant  peu  de  la  maille,  en  raison  de  leur  faible  petitesse  par 
rapport  a  celle-ci. 


Figure  4.  Mailles  amorphes  construites  et  utilisees  dans  les  calculs  de  rg.  A  gauche,  modele 
montrant  les  molecules  contenues  dans  une  unite  de  maille,  a  droite,  la  representation  dans  des 
conditions  limites  periodiques,  incluant  les  molecules  et  atomes-image. 

a)  HTPB 


d)  GAP 


d)  Gpl 


e)  Melange  GAP-Gpl 
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Ces  mailles  ont  ensuite  ete  equilibrees,  puis  soumises  a  une  hausse  de  temperature  et,  ont  fait 
Pobjet  des  calculs  de  la  dynamique  moleculaire,  sous  ensemble  NPT,  pour  diverses 
temperatures.  A  chaque  reprise,  les  50  demieres  picosecondes  de  la  simulation  servent  a 
determiner  le  volume  que  la  maille  occupera  a  la  temperature  de  la  simulation.  On  retrouve  a 
la  figure  5  les  graphiques  des  volumes  specifiques  en  fonction  du  temps  pour  les  differents 
systemes  modelises.  Dans  chaque  cas,  les  graphiques  de  correlation  de  paires 
intermoleculaires  g(r)  permettent  de  constater  la  proximite  des  divers  groupements,  puisque 
c’est  la  mesure  de  la  distribution  des  distances  interatomiques  entre  chaines  polymeres  et 
polymere-plastifiant.  Dans  les  deux  cas,  la  courbe  d’interactions  polymere-plastifiant  apparait 
a  des  valeurs  nettement  plus  elevees  que  celle  des  courbes  polymere-polymere  et  plastifiant- 
plastifiant,  indiquant  que  les  interactions  entre  deux  composantes  des  melanges  sont 
importantes,  ce  qui  est  un  indicatif  de  melanges  miscibles.  Dans  le  cas  de  la  courbe 
HTPB-DOA,  la  presence  d’un  leger  maximum  entre  6  et  8  A  peut  etre  associee  &  l’effet  d’un 
trou  de  correlation  selon  lequel,  a  certaines  distances,  les  contacts  intermoleculaires  sont 
limites  en  raison  d’un  effet  ecran  des  nuages  electroniques  qui  nuit  au  rapprochement  des 
atomes.  Le  maximum  observe  dependrait  alors  de  la  distance  de  premiers  voisins,  qui  se  situe 
aux  environs  de  5  A,  comme  on  l’a  precedemment  observe  pour  le  polybutadidne  1,4-cis  [20]. 

Figure  5:  Fonction  de  correlation  de  paires  intermoleculaires  g(r) 

a)  Melange  HTPB-DOA 
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b)  Melange  GAP-Gpl 
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3.3  Determination  des  transitions  vitreuses  des  composantes 
pures  et  des  plastifiants 


Le  but  de  notre  etude  etait  d’evaluer  si,  en  utilisant  la  modelisation  moleculaire,  on  pouvait 
faire  des  predictions  fiables  de  la  temperature  de  transition  vitreuse  de  molecules  et  melanges 
energetiques. 

Nous  avons  mene  notre  etude  en  suivant  la  methode  proposee  par  Rigby  et  Roe  [11],  tel  que 
discute  dans  la  section  2.4.  Cette  methode  est  bien  etablie  dans  les  publications  scientifiques. 
Afin  d’accelerer  les  calculs,  nous  avons  choisi  des  mailles  de  faible  dimension  :  entre  1 8,6  et 
23,4  A,  de  meme,  le  temps  de  simulation  est  reste  plutot  court.  Les  resultats  obtenus  pour  les 
composes  purs  et  les  melanges  sont  montres  a  la  figure  6.  Dans  tous  les  cas,  chaque  point 
represente  une  moyenne,  sur  50  ps,  des  valeurs  obtenues  par  dynamique  moleculaire  a  partir 
du  meme  modele.  Comme  le  nombre  d’atomes  contenus  dans  la  maille  reste  constant,  on  peut 
aisement  calculer  la  masse  volumique,  ou  sa  valeur  inverse,  le  volume  speciflque.  C’est  cette 
valeur  qui  est  trac^e  en  ordonnee  en  fonction  de  la  temperature.  Sur  chaque  graphique,  on 
peut  observer,  avant  et  apres  la  transition  a  la  phase  vitreuse,  que  le  volume  speciflque  varie 
lineairement  avec  la  temperature.  Toutefois,  les  points  sont  plus  disperses  apres  rg. 

La  difficult^  principale  consiste  a  decider  quels  points  appartiennent  a  la  droite  «  avant  T%  »  et 
a  celle  «  apres  Tg  ».  Nous  avons  utilise  la  methode  des  moindres  carres,  telle  que  mise  en 
oeuvre  dans  le  progiciel  Microsoft  Excel ,  pour  trouver  les  droites  les  mieux  ajustees  aux 
donnees  de  la  figure  6.  Avant  la  transition,  la  droite  de  regression  est  continue  et  les  donnees 
de  simulation  en  carres  pleins,  apres  la  transition  la  droite  est  en  pointilles  et  les  donnees  en 
carres  vides. 

Le  coefficient  de  correlation,  ainsi  que  la  temperature  ou  se  croisent  les  deux  droites,  qui 
correspond  Tg,  sont  enumeres  au  tableau  4.  Comme  on  peut  le  voir  dans  ce  tableau,  le 
coefficient  de  correlation  est  excellent  aux  temperatures  sous  la  transition,  dans  la  phase 
vitreuse,  (il  varie  de  0,96  a  1,00)  et  demeure  tres  bon  aux  temperatures  superieures  a  la 
transition  (de  0,88  a  0,98),  dans  la  phase  caoutchouteuse.  Puisque  la  densite  d’un  polymere 
vitreux  varie  lineairement  en  fonction  de  la  temperature,  seulement  quelques  points  sont 
necessaires  pour  d’obtenir  une  droite  avec  un  fort  coefficient  de  correlation.  En  contrepartie, 
l’etat  caoutchouteux  presente  des  variations  de  densites  plus  grandes.  Ces  variations  sont  la 
principale  source  dMncertitude  sur  la  valeur  calculee  de  rgpar  dynamique  moleculaire,  rgDM. 

A  cette  cause,  s’ajoutent  les  differences  entre  la  simulation  et  l’echantillon  reel,  notamment  sa 
teneur  en  impuretes  ou  en  humidite,  de  meme  que  les  differences  d’indices  de  polydispersite 
(ou  polymoiecularite),  /p,  qui  est  fixee  a  1,00  dans  la  simulation.  Ce  n’est  pas  le  cas  lors  de 
futilisation  de  polymeres  r£els,  sauf  s’ils  sont  obtenus  par  polymerisation  vivante.  Or,  les 
polymeres  energetiques,  ou  pour  ceux  utilises  dans  les  formulations  avec  les  molecules 
energetiques,  sont  rarement  synthetises  de  cette  maniere.  Ainsi,  le  polymere  reel  peut  contenir 
des  reticulations,  ce  qui  n’a  pas  £te  simul£  ici.  II  est  interessant  de  voir  que  fecart  entre  la 
valeur  de  T&  de  la  literature,  rglm  et  celle  obtenue  experimentalement,  Tgcxp,  qui  varie  de  3  a 
10°,  est  du  meme  ordre  de  grandeur  que  celui  obtenu  entre  la  TgSimulee,  7gDM  et  rghtt  (de  6  a 
20°).  II  est  clairement  etabli  qu’en  raison  de  la  haute  vitesse  de  refroidissement  utilisee  lors 
des  simulations  par  rapport  a  la  valeur  experimentale  (soit  10°/min),  la  temperature  de 
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transition  vitreuse  simulee  est  plus  elevee  que  celle  mesuree  experimentalement  [15].  Ceci  est 
du  au  fait  que  le  syst£me  n’a  pas  dispose  du  temps  necessaire  pour  atteindre  la  densite 
correspondant  a  un  veritable  etat  d’equilibre  thermodynamique  stable  pour  une  temperature 
donn^e.  Compte  tenu  du  temps  court  de  la  dynamique  moleculaire  de  moderation,  les  grands 
mouvements  des  chaines  moieculaires  qui  se  produisent  dans  la  realite  ne  sont  pas  observes. 
Par  consequent,  les  temperatures  de  transition  obtenues  sont  systematiquement  superieures  a 
la  veritable  Tg  du  systeme,  ce  qui  explique  l’ecart  obtenu  entre  les  valeurs  predites  et  celles 
mesurees. 

Si  on  compare  maintenant  les  ecarts  observes  entre  les  valeurs  de  TgDM  et  Tg€xp  pour  les 
composes  purs  et  les  melanges,  on  remarque  une  me  me  tendance:  l’ajout  de  DO  A  au  HTPB  a 
provoque  une  diminution  de  1 8°  de  la  Tg,  alors  que  la  dynamique  moleculaire  predit  une 
diminution  de  15°.  Des  valeurs  comparables  de  12  et  15°  sont  aussi  observees  lors  de  l’ajout 
de  GP1  au  GAP.  Ceci  implique  que,  pour  avoir  une  estimation  de  l’effet  sur  la  Tg  de  l’ajout 
d’une  substance  a  un  polymere  pur,  il  faut  modeiiser,  non  seulement  le  melange,  mais  aussi  au 
moins  une  des  substances  pures.  Dans  le  cas  present,  ces  valeurs  sont,  par  ailleurs,  de  l’ordre 
de  grandeur  de  la  precision  obtenue  par  l’estimation  du  point  de  croisement  des  droites  de 
regression  lineaire,  soit  environ  10  a  20°.  Done,  dans  l’etat  actuel  des  connaissances,  la 
modeiisation  donne  une  estimation  de  la  temperature  de  transition  vitreuse  d’une  precision 
comparable  a  celles  de  la  loi  d’additivites  ou  de  liquation  de  Fox,  surtout  lorsque  la 
difference  de  temperature  entre  les  composantes  pures  du  melange  est  faible,  comme  e’etait  le 
cas  ici.  Par  contre,  pour  des  systemes  presentant  des  interactions  interchaines  particulierement 
fortes  (interactions  ioniques,  liaisons  hydrogene),  la  modeiisation  pourrait  offrir  un  avantage 
sur  ces  methodes.  De  meme,  lorsque  les  temperatures  de  transition  vitreuses  des  composantes 
pures  ne  sont  pas  connues,  il  est  possible  de  les  modeiiser  a  partir  de  dynamique  moleculaire, 
et  il  est  alors  preferable  d’effectuer  aussi  les  calculs  des  melanges  de  la  meme  maniere. 

Tableau  4.  Comparison  entre  les  Tg  simuiees  par  dynamique  moleculaire  et  les  Tg 
experimentales  obtenues  par  nous,  et  celles  publiees. 


Syst6me 

T  «*P 

T  DM 

'g 

Coefficient  de  cortelation 

T  DM  T  exp 

/g  -/ o 

T  DM  T  litt 
fl  ”'Q 

K 

K 

Phase  vitreuse 

Phase  caoutchouteuse 

K 

K 

±10 

±5 

±  5 

HTPB 

192 

201 

1,0000 

0,9298 

9 

3 

HTPB-DOA 

177 

183 

0,9946 

0,9828 

6 

- 

GAP 

227 

240 

0,9671 

0,9665 

13 

3 

GAP-Gpl 

213 

228 

0,9620 

0,9351 

15 

10 

Gpl 

205 

225 

0,9945 

0,8842 

20 

- 
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Figure  6.  Evolution  du  volume  specifique  decoulant  des  simulations  de  dynamique 
moleculaire  a  diverses  temperatures.  Les  droites  ont  ete  obtenues  par  regression  lineaire. 


a)  HTPB  pur 


T  (K) 


b)  HTPB-DOA  (polymere  et  plastifiant) 
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(cm3  g'1) 


c)  GAP  pur 


d)  Gpl  pur 
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e)  Gap-Gpl 


T  (K) 


Recemment,  Lyulin,  Balbalev  et  Michels  [30]  ont  utilise  Pequation  de  Vogel-Fulcher- 
Tammans  pour  quantifier  l’ecart  entre  les  valeurs  de  Tg  experimentales  et  celles  determinees 
par  la  dynamique  moleculaire.  Soldera  et  Metlata  [3 1  ]  ont  propose  que  cet  ecart  suive  une 
expression  mathematique  equivalente,  la  loi  de  Williams-Landel-Ferry  ou  WLF.  Cette 
demiere  relie  le  temps  experimental  a  la  temperature  selon  : 


log  Oj  =  log  —  = 


-Ci(T-T„) 

C2+T-T0 


(7) 


Dans  F equation  (7),  aT  represente  la  variable  reduite  du  facteur  de  deplacement  necessaire  a 
la  superposition  de  deux  phenomenes  equivalents,  mais  mesures  a  deux  regimes  differents, 
tjt  est  le  rapport  entre  les  temps  de  mesure  ou  temps  de  relaxation  aux  temperatures  T  et  T<» 

T0  etant  la  temperature  de  reference  a  laquelle  les  constantes  C\  et  C2  ont  ete  etablies. 
Celles-ci  sont,  en  premiere  approximation,  constantes  pour  tous  les  polymeres,  et  des  valeurs 
de  C 1  =  17,44  et  C2  =  5 1 ,6  ont  ete  proposees  lorsque  la  temperature  de  transition  vitreuse, 
est  prise  comme  temperature  de  reference.  Ferry  indique  cependant  que  ces  valeurs 
universelles  ne  devraient  etre  utilisees  qu’en  dernier  recours,  lorsqu’aucune  valeur 
experimentale  n’est  connue  pour  le  polymere  etudie,  puisque  differents  polymeres  ont  des 
valeurs  de  C 1  et  C2  l^g^rement  differentes  a  la  Tg  [32]. 


L’ utilisation  de  ces  valeurs  standards  pour  le  systeme  actuel,  utilisant  une  valeur  de  tjt?  egale 
au  rapport  entre  les  vitesses  de  refroidissement  (en  DM)  ou  de  chauffage  (en  DSC),  conduit  au 
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calcul  d’un  facteur  correctif  (T-  Tg )  de  80°,  valeur  nettement  superieure  a  la  valeur  observee 
de  rgDM-  rgcxpdonnee  au  tableau  4,  qui  se  situe  entre  6  et  20°.  Puisque,  des  ecarts  entre  10  et 
20°  sont  aussi  notes  par  de  nombreux  autres  auteurs  [2,  4,  21],  l’utilisation  de  ^equation  WLF 
pour  corriger  les  valeurs  de  dynamique  moleculaire  n’est  pas  actuellement  conseillee. 
Plusieurs  facteurs  semblent  apporter  des  contributions  similaires  a  cette  correction.  On  a  deja 
souligne  les  differences  entre  les  polymeres  reels  et  simules  (indice  de  polymolecularite, 
tacticitd,  impuretes).  A  ceci  s’ajoute  le  fait  que  les  constantes  universelles  WLF,  C\  et  C2, 
devraient  etre  ajustees,  puisque  les  melanges  etudies  ici  comportent  de  petites  molecules 
(DOA)  et  des  oligomeres  (GP1,  GAP  et  HTPB),  pour  lesquels  les  constantes  universelles  ne 
semblent  pas  s’appliquer. 
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4.  Conclusions 


Ce  travail  avait  pour  but  de  verifier  dans  quelle  mesure  la  moderation  moleculaire 
atomistique  permettait  d'evaluer  la  temperature  de  transition  vitreuse,  Tg,  de  melanges 
energetiques.  Nous  concluons  que  Ton  peut  estimer  la  temperature  de  transition 
vitreuse  aussi  bien  des  composes  purs  que  des  melanges,  avec  une  incertitude  variant 
de  cinq  a  vingt  degres.  Cependant,  lorsque  la  temperature  de  transition  vitreuse  des 
composants  d’un  melange  est  connue  et  sa  miscibi lite  bien  ytablie,  et  que  la  difference 
de  temperatures  de  transition  vitreuse  entre  les  deux  composantes  est  faible,  la 
prediction  sur  la  base  de  la  loi  d'additivite  ou  de  F  equation  de  Fox  est  tout  aussi 
precise. 

On  peut  done  conclure  que,  compte  tenu  du  temps  de  calcul  et  la  complexity  des 
operations  effectuees  pour  arriver  a  la  valeur  de  la  Tg,  la  moderation  serait  plus  utile 
si  on  la  reserve  : 

1)  au  calcul  de  la  temperature  de  transition  vitreuse  de  substances  pures  pour 
lesquelles  on  ne  dispose  pas  de  valeurs  experimentales; 

2)  a  la  prediction  de  la  miscibilite,  tel  que  decrit  dans  Particle  [29]; 

3)  a  la  prediction  de  la  temperature  de  transition  vitreuse  lorsque  l’on  prevoit  des 
interactions  interchaines  particulierement  fortes,  comme  lorsque  des  liaisons 
hydrogene  ou  des  interactions  ioniques  sont  presentes,  puisque  les  lois  d’additivite  ou 
F  equation  de  Fox  ne  s’appliquerait  generalement  pas  dans  ces  cas  particulier; 

4)  la  prediction  de  systemes  pour  lesquels  la  difference  de  Tt  entre  les  composantes 
pures  est  grande. 

La  prediction  du  phenomene  de  transition  vitreuse  par  modelisation  moleculaire  est 
encore  en  plein  developpement.  Utilisee  judicieusement,  cette  methode  peut  s’averer 
un  outil  predictif  de  choix.  Cependant,  les  temps  de  calcul  necessaires  et  l’erreur 
associee  a  cette  technique  predictive  en  limitent  l’utilite  dans  un  contexte  industriel. 
Plus  qu’un  outil  pfedictif,  cette  technique  apporte  des  informations  permettant  de 
raffmer  notre  comprehension  des  phenomenes  mofeculaires  se  produisant  lors  de  cette 
transition.  C’est  dans  cette  optique,  que  la  methode  de  modelisation  moleculaire  est 
actuellement  presentee  dans  les  publications,  et  qu’elle  peut  apporter  des 
renseignements  pertinents  en  rapport  avec  la  transition  vitreuse. 
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Annexe  1  :  Equations  du  champ  de  force  COMPASS  (Materials  Studio) 


=  Z [Ht  (b  -  b0 )2  +  ( b -b0?  +  K< (b  - b0 )4  ] 
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Figure  2.  Equations  du  champ  de  force  COMPASS  {Materials  Studio) 

Les  termes  dans  les  equations  donnees  dans  la  figure  2  correspondent  a  : 

>  Le  terme  (1)  est  Eb.  energie  de  deformation  d’une  liaison  de  distance  interatomique  b 
par  rapport  a  un  etat  stable  d’energie  «  zero  »  definie  par  les  parametres  bG. 

>  Le  terme  (2)  est  Ee  •  energie  de  deformation  d’angle  de  valence  9 ,  est  definie 
similairement  avec  un  etat  stable  90\ 

>  Le  terme  (3)  est  EA  •  energie  de  deformation  d’angles  de  torsion  cp,  est  definie  par 

rapport  a  plusieurs  minima  possibles  (f)  i,  (j>  2  et  (f>  3.  Cette  equation  permet  de 
converger  vers  le  minimum  le  plus  proche,  par  exemple  vers  des  angles  anti  ou 
eclipses  dans  une  chaine  aliphatique,  et  permet  egalement  d’estimer  correctement  la 
difference  d’energie  entre  ces  conformations; 

>  Le  terme  (4)  est  Ez  •  terme  d’energie  hors-plan,  qui  permet  d’imposer  la  coplanearite 

de  groupements  ayant  des  electrons  delocalises,  des  cycles  aromatiques  par  exemple, 
par  le  biais  d’une  penalite  energetique,  lorsque  Tangle  x  entre  un  atome  et  le  plan  des 
autres  atomes  devant  rester  coplanaire,  devie  de  0°; 

>  Les  termes  de  (5)  a  (1 1)  sont  les  termes  croises  d’energie,  permettant  d’inclure  l’effet 
de  la  deformation  simultanee  de  deux  elements  geometriques  (ex:  Ebb%  qui 
correspond  a  la  deformation  des  distances  de  liaisons  adjacentes); 

>  Le  terme  (12)  Ec|stat:  T energie  electrostatique,  exprimee  par  la  loi  de  Coulomb,  qui 
permet  d’evaluer  l’effet  sur  Tenergie  des  charges  partielles  q{  et  q,  presentes  sur  deux 
atomes  non  lies  i  et  j  separes  par  une  distance  ry  dans  un  milieu  ayant  une  constante 
di^lectrique  z  (egale  a  1  dans  le  vide); 

>  Le  terme  (13)  E\dw:  Tenergie  non  liee  de  van  der  Waals  entre  deux  atomes  non  li£s  i  et 
j,  separes  par  une  distance  rMi 

Les  constantes  K ,  //,  F,  A  et  B ,  de  meme  que  les  valeurs  a  Tequilibre  b0t60 ,  (p° ,  <p2°et<p30,  sont 
parametrees  pour  les  differents  types  d’atomes  et  de  combinaison  d’ atomes.  L’ energie  d’une 
molecule  sera  calculee  par  la  sommation  de  toutes  ces  contributions  geometriques  entre 
atomes  lies  entre  eux  (aussi  nomme  energie  interne)  et  entre  atomes  non  lies  (ou  energie 
exteme). 
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Annexe  2:  Protocoles  cTutilisation  de  Materials  Studio  (  Version  3.2) 


1 .  Utiliser  cTabord  les  tutorials  pour  les  modules  Quick  Start ,  Visualizer  et  Amorphous 
Cell . 

2.  Dessiner  les  molecules  de  plastifiant  et  de  liant,  chaque  molecule  devant  etre  cre^e 
dans  un  fichier  .xsd  separe.  (voir  le  tutorial  de  Visualizer). 

3.  Faire  un  «  Adjust  Hydrogen  »,  et  un  «  Clean  »  pour  toutes  les  structures.  Cliquer 
ensuite  dans  le  menu  «  Build  |  Bonds  »  et  sur  «  Calculate  »  afin  que  tous  les  liens  de 
valence  soient  corrects.  Seul  des  groupements  non  conventionnels  tels  que  Fazoture  (N3) 
peuvent  presenter  des  problemes,  il  faut  dessiner  ce  dernier  de  la  maniere  suivante  : 
N=N=N. 


hfl  Bond  Calculation 


*1 

Bonding  Scheme  |  Connectivity  Options) 


j  N  o  bonding  to  s  and  f  shell 


~3 

Reload  |(aj|  t 


Options 


R  Calculate  connectivity 
R  Calculate  bond  type 
R  Convert  representation  to  |  Kekule 


r  Monitor  bonding 


Calculate 


Help 


Le  groupement  azide  est  traite  par  le  champ  de  force  COMPASS ,  mais  pas  par  PCFF. 

Les  chames  de  polymeres  doivent  etre  creees  par  le  module  Build  Polymer  (voir  le 
tutorial  Amorphous  Cell).  S’assurer  de  creer  une  chaine  avec  la  bonne  masse  moleculaire, 
correspondant  aux  donnees  experimentales  (en  ajustant  le  nombre  d’unites  de  repetition). 
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Par  exemple,  pour  creer  une  chaine  de  HTPB  avec  une  masse  de  2800  g  mol'1,  il  faut 
53  unites  de  repetition. 


fill  Homopolymer 


Polymerize  [ '.Advanced]]  Branches] 

Orientation:  [Head  to  Tail  *  1  Flip  probability.  [OO 


xJ 


Initiator:  |  hydroxyl  Terminator:  |  hydroxyl  3 

Torsion:  1 1 80  0  R7  Random 

f“  Use  random  number  seed:  [cTo 


Build 


Help 


1.  Optimisation  par  Discover  :  les  molecules  doivent  subir  une  minimisation  de  leur 
geometrie  avant  d’etre  soumises  a  la  construction  de  la  maille  amorphe  . 

a.  Choisir  d’abord  Gateway,  dans  Modules  |  Discover  |  Setup  |  Job  Control. 
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51  Discover  Setup 


Automation]  Typing] 

Energy]  Non-Bond  lob  Control 


*! 


Gateway  location:  |  logan. bioinfo. ulaval  ca_1 8888  ^  ] 

Queue:  | 

Job  description:  |  p  Automatic 

P  Run  in  parallel  onjl  — J  of  32  processors 


More... 


Help 


b.  La  minimisation  s’effectue  enfin  avec  «  Modules  |  Discover  |  Minimizer  ».  Le 
niveau  Medium  (5000  steps)  est  suffisant.  La  methode  Smart  Minimizer 
analyse  automatiquement  la  taille  du  systeme  et  selectionne  la  meilleure 
methode  de  minimisation  a  utiliser. 


1 51  Discover  Minimization 

*j 

Method: 

|  Smart  Minimizer  3] 

, 

More ... 

Convergence  level: 

(Medium  3J 

Maximum  iterations: 

[5000 

P  Optimize  cell 

Minimize  |  Files... 

Help  j 

1 
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3.  Creation  de  la  maille  amorphe  par  Amorphous  Cell. 


^Amorphous  Cell  Construction 


Construct  |  Liquid  Crystal  |  Preferences)  Setup  | 
Constituent  molecules  - 


■i 


*1 


Number 

Molecule 

Add  | 

1 

Polyc_butadiene_2793 

6 

dos 

Delete 

Temperature  (K)  Cell  type  Number  of  configurations 

j  298,0  |  Periodic  cell  1 

Density  - 

T arget  density  of  the  final  configurations:  |q  907  g/cc 

r  Ramp  density  from  an  initial  value  of:  |o"g  g/cc 

Cell  parameters  (Angstrom)  - 

a:  b:  c: 

[Specify  a,  b  H  |21 ,40  1 21  40 

W  Refine  configurations  following  construction 
r  Build  a  nematic  liquid  crystalline  phase 


Construct 


Files..  |  Help 


a)  Le  nombre  de  molecules  de  chaque  espece  doit  etre  choisi  de  fa<?on  telle  que  le  meme 
nombre  d’atomes  de  chaque  espece  soit  present  dans  le  systeme.  Ceci  permet  de  creer 
des  systemes  semblables  et  done  comparables  pour  plusieurs  paires  differentes  de 
liant  et  de  plastifiant.  Par  exemple,  a  une  chaine  de  poly  mere  de  5 14  atomes, 
correspond  6  molecules  du  plastifiant  DOA  (6  x  86  atomes  =  514  atomes).  La  density 
est  determined  avec  la  fraction  massique  de  chaque  espece,  en  connaissant  les 
densites  a  Petat  pur  de  ces  substances.  Par  exemple,  6  molecules  de  DOA 

(0,914  g/cm3)  et  une  molecule  de  HTPB  (0.,00  g/cm3)  represented  un  systeme  avec 
47,8%  de  DOA  et  52,2%  de  HTPB  (massique),  la  densite  est  done  0,478  x  0,914  + 
0,522  x  0,9  =  0,907  g/cm3. 

b)  L’onglet  Liquid  Crystal  n’est  pas  utilise  ici,  on  passe  done  aux  Preferences. 
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0  Amorphous  Cell  Construction 


Construct  [..Uguid  Crystal  ;|  Preferences]  Setup 
Mesogen  set:  JTIone" 
r  Director  vector - 

*  flji 


y:  [jxT 


,00 


“31 

Cone  angle  (deg) 

To 


xj 


Construct 


Files.. 


Help 


c)  On  ne  change  pas  les  parametres  par  defaut  de  l’onglet  Preferences  lors  des 
premiers  essais.  Si  on  rencontre  des  problemes  lors  de  la  construction  de  la  maille,  il  sera 
important  de  verifier  que  la  molecule  de  depart  a  bien  ete  correctement  minimisee;  si  ceci  ne 
regie  pas  les  problemes,  il  faudra  moduler  les  parametres  de  construction  de  la  maille. 
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fill  Amorphous  Cell  Construction 


Construct]  Liquid  Crystal  „Pre[erencesj]  Setup] 

Construction  algorithm  - 

Maximum  lookahead  bonds:  pj 

Maximum  lookahead  configurations:  pg 

Number  of  bonds  to  be  added  per  step:  Pj 


Number  of  substates  per  state:  J”j 

Substate  width:  J20 

Random  number  seed:  1 72489 

W  Check  for  ring  catenations 

Refinement  options  - 

Dynamics  steps:  pfoQQ 

Minimization  steps:  Jg 

Minimization  convergence  criterion:  |g  q-j 


*1 


Construct 


Files... 


Help 


e)  On  verifie  ici  que  le  Gateway  est  bien  «  Logan  ».  Les  deux  champs  de  force  pouvant  etre 
utilises  avec  Amorphous  Cell  sont  COMPASS  et  PCFF.  Le  champ  de  force  selectionne  ici 
determine  aussi  celui  utilise  lors  de  la  dynamique  et  lors  de  Pextraction  des  proprietes. 
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g]  Amorphous  Cell  Construction 


xj 


Construct  |  Liquid  Oystal  |  Preferences  (.'Setup]; 
Energy  evaluation 
Forcefield: 


I  compass 
[“*  Use  group  based  cutoffs 
W  Automate  energy  setup 
Job  control  - 


U 


More .. 


Gateway  location:  |  logan. bioinf o. ulaval. ca_1 8888  ^1 

- 3 


Queue: 

Job  description: 


f7  Automatic 


More...  | 


Construct  |  Files.  |  Help 


f)  La  fonction  Construction  permet  d’effectuer  une  minimisation  de  la  maille  amorphe 
(amorphous  cell)  de  «  1000  steps  »  immediatement  apres  la  construction  initiale  de  la  maille, 
il  n’est  done  pas  necessaire  de  minimiser  a  nouveau  la  maille  avec  le  Minimizer,  puisque  la 
dynamique  moleculaire  va  creer  un  etat  d’equilibre  d’energie  superieure  a  Tetat  minimise  de 
toute  fa^on.  On  cherche  ici  le  niveau  d’energie  relie  a  Tequilibre  apres  une  longue  dynamique 
moleculaire  (au  moins  100  picosecondes). 
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SJ  Amorphous  Cell  Molecular  Dynamics 


Dynamics  |  Trajectory | 

Equilibration  stage  (NVT) 
T  emperature: 


x] 


Equilibration  time; 

Production  stage 
Ensemble’ 

Pressure: 

Number  of  steps: 
Time  step- 
Dynamrcs  time: 


298,0 


1,0 


NVT 


0,0 


500000 


1,0 


500,0 


K 

ps 


il 


GPa 

i 

fs 

ps 


Run 


]  ~Bes.nl  Help 


0]  Amorphous  tell  Molecular  Dynamics 


*i 


Dynamics  I  Trajectory  l| 
Save:  [rJ 

Frame  output  every 


50000 


Run 


3 

steps 


Files.. 


Help 


Le  nombre  de  ‘frames’  (cadres)  ou  conformations  sauvegardees  dans  le  fichier  de 
trajectoire  final  est  determine  par  le  nombre  de  steps  (etapes),  divise  par  la  valeur  donnee  au 
parametre  (Frame  output  every). 


g)  La  dynamique  moleculaire  (NVT)  doit  etre  effectuee  jusqu’a  atteinte  de  l’equilibre  du 
systeme,  qui  est  represente  par  une  stabilisation  dans  la  variation  de  Penergie  potentielle  et 
l’^nergie  non  liee  ou  ‘nonbond’  du  systeme  en  fonction  du  temps. 
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Energy  vs ,  Simulation  Time 


Energy  (kal/mol) 


Legend 

Potential  Energy  —  —  Nonbond  Energy 


Systeme  non-equilibre 


Energy  vs.  Simulation  Time 


Energy  (kal/mol) 
-1200 

•1300 

-1400 

-1500 

-1600 


-1700 


-1900  - • - - - - - - - - - - - - - «— 

0  1000  2000  3000  4000 

Simulation  Time  (fs) 


Legend 

- Potential  Energy  —  —  -  Nonbond  Energy 


Systeme  a  I’equilibre 


5000 
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h)  Si  le  systeme  atteint  un  certain  equilibre,  alors  la  valeur  de  Cohesive  Energy  Density 
restera  stable  autour  d'une  seule  valeur  pour  un  systeme  donne.  Pour  calculer  cette  propriete, 
le  module  Analysis  requiert  de  selectionner  le  fichier  de  trajectoire,  puis  de  cliquer  Define 
(selectionner  la  trajectoire  desiree),  Add  to  list  et  enfin  Analyze. 


frfl  Amorphous  Cell  Analysis 


xJ 


El  Structural 
El  Energetic 

Energy  evolution 
Cohesive  energy  density 
T  emperature  profiles 
El  Fluctuation 
(S  Dynamic 
El  Mechanical 


Apply  analysis  to  - 

C  Current  structure  in 


Add  to  list 


* I 


Frame  Range 

Trajectory  Document 

ALL 

DA...  \Polyc_butadiene_2793.  xtd  [Full! 

Delete 

Clear  list  | 

T  Hip . ] 

i)  Certaines  proprietes  requierent  de  creer  des  sets,  regroupant  ensemble  des  types 
d’atomes  possedant  des  proprietes  semblables.  Pour  savoir  quel  genre  de  sets  il  faut  creer  pour 
chaque  propriete,  cliquer  le  bouton  Help  dans  le  dialogue  de  Analysis.  La  fonction  g(r)  (  pair 
correlation  function)  permet  d’analyser  1  ou  2  sets  a  la  fois,  mais  il  est  preferable  de  n’en 
soumettre  qu’un  a  la  fois  afin  de  ne  pas  se  retrouver  avec  un  trop  grand  nombre  de  courbes 
superposees  dans  un  meme  graphique. 


42 


DRDC  Valcartier  TR  2007-615 


Creation  des  ensembles  (sete) 


1 .  On  applique  d’abord  les  labels  qui  affichent  les  Forcefield  Type  de  chaque  atome. 


SI  Label 


xj 


Object  type 


Atom 


r-Font 


Custom  text 


3 


|  Helvetica  “3  p2  3 

Color  I  Ip-] 


r  Properties 


Charge 

ElementName 

ElementSymbol 

▲ 

ForcefieldTvce  ; 

FormalCharge 

Hybridization 

.j 

IsBackboneAtom 

Islsotope 

IsMixtureAtom 

1  sVisible 

zl 

Clear  Selection 


Apply 


J  Flemove  |  Remove  All  |  Help 


2.  On  selectionne  ensuite  un  type  d’atome  que  Ton  ne  retrouve  que  dans  un  seul 
compose,  par  exemple  les  carbones  de  type  c3=  des  carbones  sp2  du  HTPB. 
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3.  Ces  atomes  etant  selectionnes,  on  choisit  ensuite  I’option  Connected,  avec  Add  to  the 
existing  selection,  afin  de  s&ectionner  tous  les  atomes  \\6s  a  ce  premier  atome  s^lectionne. 


0Atom  Selection 


Select  by  Property: 


Connected 

~3 

Connected 

Bond  Count 

Attached  Bond  Type 

Forcefield  Type 

Hybridization 

Is  Backbone  Atom 

Radial  Distance 

X  Coordinate 

w 

r  Selection  mode - 

(•  Create  a  new  selection  from  all  visible  atoms 
C  Select  from  the  existing  selection 
C  Add  to  the  existing  selection 


Select 


Help 


4.  Tous  les  atomes  de  toutes  les  molecules  d’un  meme  type  etant  maintenant 
selectionnes,  on  definit  ensuite  un  nouveau  set  avec  le  Dialogue  Edit  Sets. 


New...  | 


Delete 

Select 

Help 


H  Define  New  Set 


Name  j 

W  Show  atoms 

OK  |  Cancel  |  Help  | 

L’option  Show  Atoms  n’est  pas  necessaire,  et  ne  fait  qu’encombrer  1’affichage  de  la 
maille  amorphe. 
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j)  La  connexion  au  serveur  Logan  est  necessaire  afin  de  verifier  I’etat  de  la  licence  des 
differents  modules  a  chaque  demarrage  de  Materials  Studio  (port  42333).  Les  taches  peuvent 
alors  etre  effectuees  soit  localement,  soit  sur  Logan  (port  1 8888).  Ilya  une  limite  d’  utilisation 
de  quatre  processeurs  (pour  quatre  taches  differentes)  pour  le  module 
AmorphousCellConstruct. 
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Intentionnellement  en  blanc. 
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Liste  symboles/abrev./acronymes/sigles 


MDN 

COMPASS 

PBX 

HTPB 

GAP 

DOA 

DSC 

Gpl 

NVE 

NVT 

NPT 

T 

PCFF 

CFF91 

R 


ministere  de  la  Defense  Nationale 

Condensed-phase  Optimized  Molecular  Potentials  for  Atomistic  Simulation 
Studies 

explosif  a  liant  plastique 

polybutadiene  a  terminaison  hydroxy 

Polyazoture  de  glycidyle 

adipate  de  bis(ethyl-2  hexyl) 

analyse  calorimetrique  a  compensation  de  puissance 

plastifiant  du  polyazoture  de  glycidyle 

ensemble  a  nombre  de  molecules,  volume  et  energie  constants 

ensemble  a  nombre  de  molecules,  volume  et  temperature  constants 

ensemble  a  nombre  de  molecules,  pression  et  temperature  constants 

Temperature  de  transition  vitreuse 

Polymer  Consistent  Force  Field 

Consistent  Force  Field  91 

Constante  des  gaz  parfaits. 
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